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Recently, many new nature-inspired optimization algorithms have been 
introduced to further enhance the computational intelligence optimization 
algorithms. Among them, lightning search algorithm (LSA) is a recent 
heuristic optimization method for resolving continuous problems. It mimics 
the natural phenomenon of lightning to find out the global optimal solution 
around the search space. In this paper, a suitable technique to formulate 
binary version of lightning search algorithm (BLSA) is presented. 
Three common probability transfer functions, namely, logistic sigmoid, 
tangent hyperbolic sigmoid and quantum bit rotating gate are investigated to 
be utilized in the original LSA. The performances of three transfer functions 
based BLSA is evaluated using various standard functions with different 
features and the results are compared with other four famous heuristic 
optimization techniques. The comparative study clearly reveals that tangent 
hyperbolic transfer function is the most suitable function that can be utilized 
in the binary version of LSA. 
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1. INTRODUCTION 

In real-world condition, most of the problems are highly nonlinear. Over the last few decades, 
nature-based optimization algorithms have become powerful optimization tools for solving high dimensional 
nonlinear real-world problems in the area of science, commerce and engineering [1-3]. As no optimization 
algorithms can guarantee best result for all kinds of practical problems [4] the mimicking sense from nature 
always inspire to develop new optimization technique. One of the most recent nature-inspired optimization 
algorithms is the LSA first presented by Shareef et al. [5] based on natural occurrence of lightning. 
The basic LSA chiefly works in continuous search environment to resolve real area constrained 
optimization difficulties. Nevertheless, numerous optimization problems are vital to set in binary manner and 
solve them in binary search space by choosing “0” or “1” as well as seldom the outcomes need to be 
exhibited with some bits. Therefore, binary optimization algorithms are needed to solve those problems. 
In practical, the binary optimization algorithms are successfully utilized in many applications such as feature 
selection and modeling [6], unit commitment problems [7], power quality monitor placement [8], micro-grid 
operation scheduling [9], power quality and reliability enhancement [10], flow shop scheduling [11], 
optimum power flow problem [12], knapsack problem [13], optimal placement of sensors [14] and so on. The 
typical continuous domain optimization algorithm can be converted to binary one. The key factor between 
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the continuous and binary domain optimization algorithms is the transfer function. This function maps the 
binary search space from continuous one depending on the function’s probability condition. The most 
common transfer functions, namely, log sigmoid, tangent hyperbolic sigmoid and quantum bit rotating gate 
are generally applied for updating process of the continuous version of algorithms. The log sigmoid 
activation function has been successfully applied in many binary domain optimization algorithms such as 
binary particle swarm optimization (BPSO) [15], binary firefly (BFA) [11], binary cuckoo search [16] and 
binary bacterial foraging [17]. Similarly, tangent hyperbolic sigmoid activation function is utilized to form 
binary gravitational search algorithm (BGSA) [18], binary artificial bee colony [19] and binary 
magnetic optimization [20]. Nonetheless, quantum bit is promoted on existing binary algorithm for 
performance improvement such as quantum-inspired binary particle swarm optimization [7], 
quantum-inspired binary gravitational search algorithm [21], quantum-inspired binary firefly algorithm [10], 
quantum-inspired tabu search algorithm [22] and etc. The aim of this study is to formulate a binary version of 
LSA and identify which transfer function is most suitable to modify traditional LSA to handle binary 
optimization problems. Finally, the performances of those transfer functions are distinguished in terms of 
correctness of results and conjunction speed. 


2. BINARY LIGHTNING ALGORITHM (BLSA) 
The binary lightning search algorithm (BLSA) is first introduced by Islam et al. [23] to elevate the 
binary combinational difficulties relaying on the binary search space. In BLSA, three projectile forms are 


used: the transition projectiles that generate the first step leader population, the space projectiles ( re ) that try 


to become the leader, and the lead projectile (p") that signifies the projectile excited from the best-positioned 
step leader. The updating process of space projectile is presented as 


Preis = Pp: + exprand(u, ) (1) 


where exprand is an exponential random integer, and y; signifies the distance among the lead projectile p” 
and the space projectile p under study. 


The new position of p’ at step+/ can be expressed as 
eee oe d (2) 
Pi. new = P; +normrand (Lt, Oo, ) 


where normrand is a random number caused by the normal distribution function, fz is occupied as current 
value of p’, and ois the scale parameter that exponentially declines as it discover the best solution. The full 


BLSA can be found in [23]. The following subsections, explains their adaptation in making the binary 
version of LSA. 


2.1. Logistic sigmoid assisted updating 

In this paper the probability transfers function T;(p;), taken as a probability gate, is calculated based 
on (3) to plan a binary search space. Conversely, the projectile’s position at step+/ is upgraded succeeding 
the probability function with a terms as shown in (4). This formulation must exhibit less possibility 
to alter the projectile position if | p, |1s small. It also should be capable to make great possibility of altering 


the projectile location with bigger value of| p,|. Logistic sigmoid transfer function has this feature 


as shown in (3). 


1 (3) 
T,(p;) = ——— 
(Pd lte? 
_Jpi. & rand <|T, (p,)| (4) 
vo P;, otherwise 


where rand is the uniform arbitrary variable in interval [0, 1]. 


2.2. Tangent hyperbolic sigmoid assisted updating 
Like the logistic sigmoid probability transfer function, Tangent hyperbolic activation function can 
be utilized in the updating process as shown in (5) to drawn to a binary search universe. Conversely, the 
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projectile’s position at step+/ is upgraded considering the probability function with a terms as shown in (6). 
In (5) suggests that when the | p, | are large from a specific position, the probability of directing p; in (6) to a 


new location at step+/ is high, whereas by decreasing the | p, | at later iterations, the probability of directing 


piis reduced and lastly when the | p, | 1s zero, the position of p; remains unaffected. 


T,(p,) =| tanh (p,)| ©) 
= Dis if rand <| T, (Pi) | (6) 
oe p;, otherwise 


where rand is the uniform random variable in interval [0, 1]. 


2.3. Quantum bit assisted updating 

A quantum bit (Q-bit) is the lowest unit in quantum computing in which it stores the data in dual- 
state. This Q-bit can be in the “0” state, in the “1” state or in a linear superposition of two [24]. The state of 
Q-bit can be represented as follows: 


lv) =a|0) + HII) ) 


where @and £ are complex numbers that fulfill the probability magnitudes of the resulting states. Hence, 
the states can be normalized to 1 as shown [24]: 


leP +a) (8) 


The Q-bit state is upgraded by using a quantum gate, known as alterable gate, and can be characterized as a 
unitary operator U. 


cos(A@,) —sinCA@,) | (9) 


7 a be A@)  cos(A@) | 


where A@,, i= 1,2,3,...,n denotes the rotation angle of each Q-bit toward either the “0” or “1” state depending 


on the Q-bit sign. The rotation angle can be expressed as 
KG ps, (10) 


where @is the rotation angle amplitude beside the iteration and decline monotonously from O,,,, 10 O.,,, a8 

shown in (11). The Q-bit individual sequence is then upgraded based on the rotation angle and 
2 

rotation gate (12). Finally, the projectile’s position is upgraded based on the probability of Z| 


as shown in (13). 


= npg — ypare — Oppin) ee Pesrtat (MH) 
STEP max 
@(step +I) | = 50. 6) Pinel (12) 
|B (step + 1 | | B,Cstep) 
a _ fl, & rand <|£B, (step +) (13) 
v_new “" 10, otherwise 


where rand is the even arbitrary variable in interval [0, 1]. 
The following are the steps of BLSA algorithm. 
Step 1: The algorithm is initialized by declaring the iteration, channel period and leader tips energies. 
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Step 2: Arbitrarily create step leaders (transition projectile). 

Step 3: Evaluate the performance (projectiles energies). 

Step 4: Update leader tip energies as well as the best and worst step leaders. 

Step 5: Check maximum channel time. If maximum channel time occurs, replace the worst step leader with 
the best one. 

Step 6: Then, update each projectile position using (4) or 6 or 13 and evaluate the performance. 

Step 7: Check focking occurs or not after that update step leaders. 

Step 8: Repeat Steps 4 to 7 until the optimization criteria is achieved. In this study, convergence criterion is 
based on maximum iteration number. 

Step 9: Print solution. 


3. PERFORMANCE EVALUATION USING BENCHMARK FUNCTIONS 
3.1. Benchmark function 

Benchmark function plays an important role for testing and validating new optimization algorithm. 
It governs the competence and effectiveness of the algorithm. To measure the effectiveness of BLSA, a range 
of benchmark functions are efficiently applied. It contains different search spaces, dimensions and global 
optimal are shown in Tables | to 3. The benchmark functions are classified as follows: 


3.1.1. Unimodal functions 

The main characteristics of benchmark functions F1 and F2 used in this study are high dimensional 
functions and contain only one global optimal for each function. Table 1 represents the properties of utilized 
unimodal functions in detail. 


Table 1. Unimodal test functions 


ID _ Test function Name of Function Search range _Dimension(n) _ Best known 
nh 
= 2 
no SOO=D 0x; Sphere [-100,1007 39 0 
i=l 
n n 
pp FOO=Dibel+ TT schweter2.22 [-10,10 ay 
i=l i=l 


3.1.2. High dimensional multimodal functions 

The high dimensional functions F3 and F4 are considered as the most challenging problems 
because it contains many local minima. Table 2 presents the properties of applied high dimensional 
multimodal functions. 


Table 2. High dimensional multimodal test functions 
ID _ Test function Name of Function Search range Dimension — Best known 


_ n 2 
RB S&)= Dk —10cos(27x,) + 10] Rastrigin -5.12,5.12" 39 0 


FCx) 20e%0{ 0.2 | ey 
nN j=A 


exo LS cos(27x; >} +20+e 
Nj 


30 0 


F4 Ackley [-32,32]" 


3.1.3. Binary functions 
This function group is basically consisted of binary domain test function where global optimum 
depends on the function dimension. Table 3 reveals the properties of utilized binary test functions. 
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Table 3. Binary test functions 


ID Test function Name of Function Search range Dimension (n) _ Best known 
pp S(O) = 2 oF Max one [op 80 80 

7 8i 
6 FCO) =D>IC TI] x) Royal road [op = 

i=l) j=8CG—1+1 


4. PERFORMANCE TESTS 

In this section, the performance of proposed BLSA with various bit updating procedures is shown 
using different benchmark functions and compared with other four heuristic optimization techniques to 
validate the results. 


4.1. Parameter settings 

Like other heuristic optimization techniques, BLSA has its own algorithm dependent parameter for 
its proper functioning. The performance of BLSA is generally influenced by the energy Ept. In this study, the 
value of Ept is set as 1.05. To validate the results, the effectiveness of BLSA is validated with BPSO [15], 
BFA [11], genetic algorithm (GA) [25] and BGSA [18]. To make fair comparison, the population size and 
maximum iteration are set as 50 and 500, respectively, for all algorithms. 


4.2. Test results and discussion 

The simulations are conducted to solve the above-mentioned benchmark functions under various test 
function groups in order to evaluate the performance of BLSA with each transfer function in which BLSAlog 
(BLSA_L), BLSAtan (BLSA_T) and BLSAQbit (BLSA_Q) represent the log sigmoid, tan hyperbolic 
sigmoid and Quantum bit, respectively, based BLSA. Among three test cases, Tests 1 and 2 are related to 
objective function minimization whereas Test 3 belongs to objective function maximization. 


4.2.1. Test 1 

This test group is consisted of unimodal benchmark functions F1 and F2 which are high dimensional 
in nature. The functions Fl and F2 are comparatively easy to solve because of single global optimal. Each 
benchmark function is executed 30 times independently for all algorithms while the dimension of these test 
functions is set as 30. A comparative study is shown in Table 4 in terms of best, worst, average, median, 
standard deviation and minimum number of iteration. The best performance for each algorithm is weighted in 
bold face. From Table 4, it is shown that the BLSA can find the best known value with all bit updating 
mechanisms. However, BPSO, BFA, GA and BGSA cannot find the global optimal solution. The results of 
these algorithms are further illustrated in box plots as shown in Figure 1 relied on 30 independent runs and 
the convergence characteristics are depicted in Figure 2. From Table 4 and Figure 2, it is clearly observed 
that the tangent hyperbolic sigmoid activation function based BLSA (BLSAtan) required less number of 
iterations to reach the global optimal solution. In some cases, BPSO as shown in Table 4 requires less 
number of iterations, however, it cannot find the best known value. 


Table 4. Global optimization results of Test 1 for benchmark function in Table 1 


Test Fn. — Algorithm Best Worst Average Median Standard Deviation Minimum iteration 
BLSAtog 0 0 0 0 0 62 
BLSAtan 0 0 0 0 0 28 
BLS Aapit 0 0 0 0 0 34 

Fl BPSO 20239.97 49226.20 33774.02 34494.64 7030.68 23 
BFA 34966.2 63715.34 52878.02 54068.65 7665.34 297 
GA 46584.69 75625.58 66088.58 66521.42 6120.112 364 
BGSA 1.53E-05 17.03518 1.196955 0.012367 3.292828 373 
BLS Aiog 0 0 0 0 0 94 
BLSAtan 0 0 0 0 0 20 
BLS Aanit 0 0 0 0 0 36 

F2 BPSO 50.078 19 102.6368 79.53768 78.85152  13.00522 21 
BFA 84.82404 139.8556 117.5649 118.5391 — 11.91653 371 
GA 129.2736 161.879 146.3888 147.5646 = 10.14098 56 
BGSA 0.011719 4.082031 0.598079 0.259764 0.874547 385 
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Figure 1. Global optimization results of Test 1 for benchmark functions F1 to F3. (a) Fl and (b) F2 


x10" FA 


a 
va 


Function value 
o 


—— oe oT 


i-- 
-— 


0 100 200 300 
Iteration 


(a) 


400 500 


Function value 


0 100 


200 


Iteration 


(b) 


300 400 500 


Figure 2. Convergence characteristic curves of Test 1. (a) Fl and (b) F2 


4.2.2. Test 2 


In this test, exploration and exploitation abilities of the BLSA are tested using multimodal high 
dimensional functions. It is important to note that the complexity level of the problems have increased 
because of containing many local minima. Due to several local minima, it is usually hard to find the best 
optimal solution. Therefore, final optimal outcomes of the tested algorithm are more significant. Similar to 
Test 1, again each benchmark function is performed 30 times while function dimension is maintained as 30. 
The numerical results are tabulated in Table 5. The best performance for each algorithm is weighted in 
boldface. The outcomes of 30 independent runs are depicted in Figure 3 through box plots. Furthermore, 
the convergence characteristic curves are shown in Figure 4. Similar to previous test, the tangent hyperbolic 
sigmoid activation function updating procedure (BLSAtan) is faster to find the global optimal solution than 


the others. 


Table 5. Global optimization results of Test 2 for benchmark function in Table 2 


Test Fn. Algorithm Best Worst Average Median Standard Deviation | Minimum iteration 
BLSAtog 0 0 0 0 0 74 
BLSAtan 0 0 0 0 0 21 
BLS Aanit 0 0 0 0 0 30 

F3 BPSO 232.9097 403.0246 320.4552 323.8423 42.55152 28 
BFA 278.638 444.8202 387.1754 392.0954  39.73813 477 
GA 359.775 474.774 436.5995 438.7085 25.06685 44 
BGSA 3.811253 15.94872 8.651355 8.373049 +=. 2.764086 384 
BLSAtog 8.88E-16 8.88E-16 8.88E-16 8.88E-16  1E-31 52 
BLSAtan 8.88E-16 8.88E-16 8.88E-16 8.88E-16  1E-31 28 
BLS A anit 8.88E-16 8.88E-16 8.88E-16 8.88E-16  1E-31 34 

F4 BPSO 19.56376 20.72198 20.34145 20.44813 0.310445 17 
BFA 20.34612  20.72151 20.54059 =20.55727 (0.105571 227 
GA 20.6778 21.08882  20.96458  20.96271 0.08477 2 
BGSA 0.144861 —11.39533 3.744453 2.980169 2.872996 385 
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Figure 3. Global optimization results of Test 2 for benchmark functions F3 and F4. (a) F3 and (b) F4 


600 : : 
---, ---BLSA,, 
500 ===BLSA, 
7 ===BLSA,,, 
I === BPSO 
400k | itqesm iia cee th 
2 i ee ee bbe 4 
$ tot x =="GA A 
£3007! I 1 ===BGSA [I + 
5 . I a =a 
5 ee Se ee eee eee Sen 
Z rv ~ 
20011 | I .4 
I I \ 
I I V 
1004 I 1 *. 
| 
| 
* 


0 100 


200 


Iteration 


(a) 


400 


500 


ial =“ TTT Teeeeeeeeeee 
soll \ i ueneneeieie | i iii 3 
i! \ 
® il ay 
2 15/yi | \ 
5 jul \ ==-BLSA 
< wil 1 log 
fist * ==-BLSA, 
e tory! ! 4 BLSA,,,| | 
2 jill 1 ——-BLSA,, 
J H 1 === BPSO 
stl I 1 === BFFA 
il 1 \ ===GA 
I ===BGSA 
il I \ 
a 
0 100 200 300 400 500 
Iteration 
(a) 


Figure 4. Convergence characteristic curves of Test 2. (a) F3 and (b) F4 


4.2.3. Test 3 


To determine the capability of the BLSA for maximization problem, this test group is handled with 
binary nature (Max one and Royal road) test function. The dimensions of these functions are set as 80. 
The simulation results are summarized in Table 6 based on 30 independent runs. The best performance for 
each algorithm is weighted in boldface. Similar to Tests 1 and 2, the performance of BLSA are identical for 
various bit updating procedures with different transfer functions as shown in Table 6 and Figures 5 and 6. 
For maximization problems, convergence characteristic curves portrayed in Figure 6 clearly represents that 
the quantum bit assisted BLSA (BLS AQbit) outperforms the others. 


Table 6. Global optimization results of Test 3 for benchmark function in Table 3 


Test Fn. Algorithm Best Worst Average | Median Standard Deviation Minimum iteration 
BLSAtog 80 80 80 80 0 104 
BLSAtan 80 80 80 80 0 36 
F5 BLSAapit 80 80 80 80 0 13 
n=80 BPSO 71 65 67.6 68 1.499425 30 
BFA 74 70 71.93333 72 1.048261 227 
GA 70 65 66.73333 66 1.201532 320 
BGSA 80 80 80 80 0 302 
BLSAtog 10 10 10 10 0 110 
BLSAtan 10 10 10 10 0 41 
F6 BLSAauit 10 10 10 10 0 15 
n=80 BPSO 5 2 3.466667 3 0.62881 23 
BFA 6 4 4.233333 4 0.504007 69 
GA 5 3 3.5, 3 0.572351 197 
BGSA 10 9 9.566667 10 0.504007 295 
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Figure 5. Global optimization results of Test 3 for benchmark functions F5 and F6. (a) F5 and (b) F6 
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Figure 6. Convergence characteristic curves of Test 3. (a) F5 and (b) F6 


5. CONCLUSION 

This paper investigated different methods to realize binary version of LSA using logistic, tangent 
hyperbolic and quantum bit probability gate, respectively. These activation functions modified the position 
update procedure of the original version of LSA. The solution finding ability and consistency of various 
transfer functions assisted BLSA were tested with various benchmark functions of different complexity. 
The performance of BLSA is compared with other heuristic optimization techniques. The comparative study 
reveals that the performance of BLSA with tangent hyperbolic sigmoid function projectile and step leader 
updating procedure is the best to be utilized for minimization problems whereas quantum bit method is better 
for maximization problems. Moreover, it can also be concluded that both tangent hyperbolic and quantum bit 
assisted BLSA can be equally used for both minimization and maximization problems since they can provide 
best global optimal solution for all standard benchmark functions corresponding to top known value without 
significant differences in convergence speed and computational burdens. 
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